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Abstract 

Temperature-independent magnetic viscosity in ferritin has been observed 
from ~ 2 K down to 100 mK, proving that quantum tunneling plays the main 
role in these particles at low temperature. Magnetic relaxation has also been 
studied using the Landau-Zener method making the system crossing zero res- 
onant field at different rates, a = dH/dt, ranging from 10 -5 to 10 -3 T/s, and 
at different temperatures, from 150 mK up to the blocking temperature. We 
propose a new Tln(AH e tf /tqol) scaling law for the Landau-Zener probability 
in a system distributed in volumes, where AH e ff is the effective width of the 
zero field resonance. 
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Over the past decade there have been experimentally observed a large number of quantum 
phenomena in the dynamics of the magnetic moment of mesoscopic systems. Monodomain 
magnetic particles of nanometer size contain thousands of magnetic atoms strongly intercon- 
nected by exchange interaction. As a result of the exchange interaction, all the atomic spins 
align parallel or antiparallel between them, resulting in a ferro, or ferri-antiferromagnetic 
ordering respectively. One particular system has had much attention since Awschalom and 
co-workers announced the observation of a resonance near 1 MHz that interpreted in terms 
of quantum coherence of the magnetic moment [|l]||: the system is composed by antiferro- 
magnetic particles which grow inside the cage of the horse spleen ferritin proteins 0. 

Next experimental studies of the dynamics of the magnetization of ferritin particles, 
carried out at kelvin regime, showed different phenomena interpreted as quantum tunneling 
of the ferritin magnetic moment |3|-[T0|. These phenomena can be differentiated as follows: 
The first was the temperature-independent relaxation of the magnetization below 2.3 K 
observed in the magnetic viscosity measurements [Q. Secondly, a non-monotonic behaviour 
of the blocking temperature, T B , on the magnetic field B,|^,^|-|IT| , contrasting with the 
monotonic classical expectation. Thirdly, a maximum at zero field in the magnetic field 
derivative of the magnetization extracted from the hysteresis cycles over 4 K 0,|5|-Q. And 
finally, a rapid increase of the magnetic viscosity as the magnetic field approaches to zero at 
temperatures higher than 3 K |||8|||. More recently, there have been done 57 Fe Mossbauer 
spectroscopy measurements on an artificial ferritin sample down to 50 mK showing non- 
incoherent tunnel fluctuations around 10 8 Hz ||1 1|| . Attending to the recent observation 
of resonant quantum tunneling of the spin observed in molecular clusters like Mn 12 -Acetate 



i4|T3|l , one can conclude that all these phenomena observed in ferritin particles in the kelvin 
range can be attributed to thermally activated resonant quantum tunneling of the magnetic 
moment at zero field. The observation of tunneling only at zero field lies in the fact that 
ferritin has a broad distribution of energy barriers due to the distribution of volumes of 
the particles (f(U = KV), where K is the anisotropy constant), and their anisotropy axes 
are randomly oriented. In the other hand, the temperature-independent viscosity below the 



crossover temperature, T c = 2.3 K, indicates that, below this temperature, the quantum 
magnetic relaxation of ferritin occurs through the lowest states. 

In this paper we present new magnetic data which extend the quantum relaxation mea- 
surements to the millikelvin regime. At the same time, in order to estimate the value of the 
quantum splittings of ferritin particles, we have done measurements of the change of magne- 
tization when the system crosses zero magnetic field at different rates of the field sweep and 
we analyze the results in terms of the Landau-Zener probability associated to the magnitude 
of the splitting playing the main role in the quantum relaxation. This is the same method 
used by Wernsdorfer et al. to determine the quantum splittings of Fes molecular clusters 
0. 

Ferritin is an iron storage protein. It has a spherical cage of about 8 nm in diameter in 
whose interior grows mineral ferrihydrite combined with a phosphate. Its core is equivalent 
to a small antiferromagnetic particle. The size of the core in natural ferritin ranges from 3 to 
7.5 nm. The fully packed ferritin contains 4500 Fe 3+ ions. A small magnetic moment of the 
particle arises from the non-compensation of collinear spin sublattices due to the finite size 
and irregular shape of the core. The spin of the sublattice, S, is of the order of 5000, while 
the non-compensated spin, s, is below 100. This number corresponds to 15 non-compensated 
Fe 3+ ions extracted from magnetic susceptibility measurements [§J, in agreement with the 
theoretical expectation for the square root of the number of ions at the surface of the 
particle, Ng ~ (4500) 2 / 3 , which gives (4500) 1 / 3 ~ 16. This non-compensated spin looks in 
one of two directions along the anisotropy axis of the particle. In our experiments we have 
used a Fluka Biochemical diluted natural ferritin sample. The distribution of volumes of 
the sample, f(V)dV, is plotted in the inset of figure 1 (extracted from reference ||). The 
center volume is Vq ~ 150 nm 3 . The low moment of the antiferromagnetic particles makes 
the interactions between different particles to be negligible, as one can see in the inverse 
susceptibility in the superparamagnetic regime that extrapolates to zero at T^0 0. 

Low temperature magnetic relaxation measurements were done in an Oxford Instruments 
3 He- 4 He dilution cryostat in the following manner: The sample is cooled until the measure 
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temperature and then a magnetic field of 1 T is applied during 10 minutes. After that, the 
field is switched off and the magnetization is measured during 3 hours. In order to avoid 
remanent fields in the superconducting magnet and to obtain the relaxation measurements as 
close as possible to zero field a demagnetizing cycle is immediately applied after switching off 
the field. The demagnetizing cycle was previously tested in a pure Pb diamagnetic sample. 
This method makes the field along the major part of the relaxation to be zero with a precision 
of ±1 Oe. The measurements were done at different temperatures ranging from 100 mK up 
to 1,5 K in the dilution cryostat and repeated in the same manner in a Quantum Design 
MPMS magnetometer at temperatures up to 25 K. The logarithmic on time dependence 
of the magnetic relaxation is clearly observed over the whole measure. In a sample with a 
distribution of energy barriers, the quantitative magnitude which measures the relaxation 



time is the magnetic viscosity defined as QT5 



s= 1 dM(H,T,t) 

(M ini (H, T) - M eq (H, T)) dln(t) ' 1 ' 

where M eq (H,T) is the equilbrium magnetization of the system at fixed temperature and 
field, which is M eq (H,T) = in our case, and Mi n i(H,T) is the initial magnetization. 
In our experiments M ini (H,T) was taken from the extrapolation at small time of each 
magnetic relaxation curve. It is known that after switching off the field the system rapidly 
runs to a critical state in a time much more shorter than the times involved in the slow 
relaxation process occurring after the system reaches this critical state and relaxes to the final 



equilibrium state [15 . The observed dependence of the magnetic viscosity with temperature 
is shown in figure 1. The viscosity shows a maximum at Tg ~ 10 K. This is the blocking 
temperature, defined as 

ln(t m /r ) 

which for viscosity measurements, with a characteristic measuring time, t m , of hours, corre- 
sponds to the unfreezing of the magnetic moment of a particle of volume V which changes 
its orientation jumping over the energy barrier. In a sample distributed in size (see inset 



of figure 1) there is a distribution of energy barriers, f(U = KV). The rate at which indi- 
vidual moments of the particles jump across the anisotropy barrier depends on temperature 
through the Arrhenius exponential factor, exp^—KV/KsT). The maximum observed in the 
viscosity at Tg corresponds to the unfreezing of the particles having the center volume, Vq, 
of the volume distribution. If we look at the volume distribution of figure 1, the blocking 
temperature may correspond to the particles with a volume around 150 nm 3 . For these 
particles, using eq. (2), with r ~ 10 -8 s, t m = 10 4 s and K = 2.5 x 1CT 5 erg/cm 3 , we 
obtain T# ~ 10 K, in good agreement with the experimental result. As the temperature 
decreases, the magnetic viscosity goes to zero, as expected for thermal relaxation in a sys- 
tem with barriers distribution. However, below ~2K the viscosity becomes independent on 
temperature down to 100 mK. This temperature, at which the system crosses from thermal 
to quantum relaxation regime is called crossover temperature, T c Q. The new data showed 
in this paper extend the observation of the plateau of the magnetic viscosity down to a few 
millikelvin. This takes high relevance assuming the fact that below T c the system relaxes 
exclusively through the lowest levels of the magnetic structure by quantum tunneling. This 
temperature does not depend on the volume of the particles. The expression expected from 
theory which determines this temperature for antiferromagnetic monodomain particles is 
T c ~ (2e an e ex ) 1 ^ 2 /2tt []16| . Taking e an ~ 0.1 K (anisotropy energy per spin) and e ex ~ 10 3 K 
(exchange energy per atom) ||, we obtain T c ~ 2 K in good agreement with the experimental 
value. 

ZFC magnetization measurements have been done from 100 mK up to 25 K. The mea- 
surements were done as follows: First the system is cooled from room temperature down to 
100 mK in the absence of magnetic field. Due to the random orientation of the anisotropy 
axes of the particles the total magnetization is zero. Then a small magnetic field, H = 
200 Oe, is applied and the magnetization is measured as the temperature is increased. The 
result is shown in figure 2. In this figure we can observe the maximum corresponding to the 
blocking temperature of the sample, T B ~ 13 K, in agreement with eq. (2) using t m = 10 
s. A significant fact that one can observe in this figure is the increase of the magnetization 
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as the temperature goes to zero below ~ 0.5 K. The behaviour of the magnetization in this 
temperature range is 1/T, indicating that there is a superparamagnetic contribution to the 
total magnetization. To understand this let us see the expected behaviour for a sample 
with a broad distribution of barriers in the thermal regime. After the field is applied at 
the lowest temperature, the smallest particles, having the smallest barriers, become free to 
rotate their magnetic moments out of the anisotropy axis and align them with respect to 
the applied field, contributing to the net magnetization of the sample. As the temperature 
increases, bigger and bigger particles can jump across the anisotropy barrier and the total 
magnetization increases. When the major part of the magnetic moments of the particles 
are free, over Tg, the system behaves paramagnetically following the 1/T Curie behaviour. 
This is opposed to the observed decrease of M(T) at low temperature. Two different facts 
can explain this paradox: (a) thermal superparamagnetic behaviour of a second distribution 
of particles of smaller size or (b) quantum superparamagnetic behaviour of the smallest 
particles of the mono-distributed sample. The first explanation does not agree with the 
distribution showed in figure 1, where there is not a significant number of particles below 50 
nm 3 . Applying eq. (2) we obtain that a significant number of particles smaller than 5 nm 3 is 
needed to obtain Tg < 0.1 K. The ac-susceptibility measurements (window time of 10~ 3 s) 
|| and Mossbauer espectroscopy (10~ 8 s) [0,0 show that there are not significant particles 
of this size behaving superparamagnetically at low temperature. The second explanation, 
quantum superparamagnetism, explains better the 1/T behaviour at low temperature. This 
is the quantum behaviour of the particles for which Tg is smaller than T c . That is, these 
particles do not feel the anisotropy barrier because they can rotate their magnetic moments 
by quantum tunneling even if the temperature is not enough to jump across the barrier. 
The quantum tunneling rate is determined by the WKB exponent, B ~ KV/ KbT c , in the 
following manner, T ~ exp(-B). This means that the particles having smallest size have 
the higher probability to tunnel across the barrier. 

The most direct way to measure the quantum tunneling splitting, A, is by using the 
Landau-Zener model [0, which gives the tunnel probability, P, when a resonance is crossed 



at a given sweeping rate, a: 

hA 2 



P = 1 — exp 



(3) 



2g\i B Sa_ 

where h is the Plank's constant and S is the spin of a particle. Due to the distribution of 
volumes in ferritin there are a distribution of spin values, S(V), and a distribution of quan- 
tum splittings, A(V). Also, the random orientation of the anisotropy axis of the particles 
in the sample introduces a distribution of sweeping rates, a(9), on the angle between the 
applied field and the anisotropy axis of each particle. This makes that different particles 
have different tunnel probability at a given sweeping rate depending in both volume and 
orientation respect to the applied magnetic field. Taking into account the mentioned condi- 
tions, we can express the change of magnetization of the whole sample as the zero resonant 
field is crossed from Hi to Hj at a given a in terms of the Landau-Zener probability as 
follows: 

Mf-Mt _ W 2 



M eq - Mi 



2tt sin{9)d9 [ S[V]P[A{V),S{V),a{9)}f[V}dV , (4) 
Jo Jv 



where Mj, Mf and M eq are the initial, final and equilibrium magnetizations, respectively. 
The integral over 9 has been chosen to take into account the random orientation of the 
anisotropy axes of the particles respect to the applied field. The form of a(9) for one 
particle is then cos{9)a. 

Our experiments where done in the following manner: First, a saturating magnetic field 
was applied at the measure temperature. Then, the field was changed to Hi = 250 Oe at the 
highest sweeping rate and the magnetization was measured giving Mj. Immediately, the field 
was changed to H f = -250 Oe at a given a, measuring Mf after the process was finished. The 
procedure was repeated at different sweeping rates, ranging from 10" 5 T/s up to 10 -3 T/s 
and at different temperatures, from 100 mK up to the blocking temperature. The results are 
shown in figure 3. In order to make the nomenclature shorter we will use Pam (probability 
to change the magnetization) instead the expression given in eq. (4). One can see that, at 
a given temperature, Pam increases when a decreases. That is, as the zero field resonance 



is crossed slower the probability to change the magnetization of the sample is higher. With 
the same dependence in a, the probability becomes higher for higher temperatures. The 
bahaviour of Pam on 1/a is perfectly logarithmic. This dependence reminds the behaviour 
of the time magnetic relaxation observed in this sample and, in general, in any sample 
with barriers distribution. Indeed, we can find the equivalence between the sweeping rate 
and time using t = AH /a, where AH = Hi — Hf = 500 Oe in our experiment. Due to 
this equivalence we can define a new parameter, Slz, to evaluate the characteristics of the 
magnetic relaxation in a Landau-Zener process with a barrier distributed sample, in the 
same manner that the magnetic viscosity does it in time magnetic relaxations. That is, Slz 
can be expressed as 

_ dP A M /j,x 
bhZ ~ d\n(AH/a) ' (5) 

where Pam = (Mf — Mi)/(M eq — Mi). The temperature dependence of the Landau-Zener 
viscosity, Slz, is shown in the inset of figure 4. From the comparison of this result with 
the magnetic viscosity extracted from time magnetic relaxations (figure 1), the agreement 
between the results of both methods is clearly observed. S L z has a maximum at 10 K, which, 
using eq. (2), corresponds to the blocking temperature with an effective time, AH /a, of 
10 4 s, this is, a = 10~ 5 T/s and AH = 5xl0~ 2 T. That is, the Landau-Zener procedure 
carried out in a sample of particles distributed in volume gives the same information that the 
magnetic viscosity analysis. However, we can extract new information from this procedure 
if we analyze the change of magnetization as the zero field resonance is crossed under a new 
scaling law proposed following. 

To look for evidences of thermal or quantum relaxation from time magnetic relaxation 
experiments it is usually used a Tln(t/r ) plot. In the thermal relaxation regime the de- 
pendence of the magnetization on Tln(t/r ) scales in a master curve if the characteristic 
relaxation time, r , is adequately chosen. This analysis permits to extract evidences of 
quantum tunneling when the magnetic relaxation departs from the master curve below a 
characteristic temperature. This temperature corresponds to the crossover temperature, T c , 



explained before. Independently, it is possible to find the barriers distribution function from 
the derivative of the master curve. We show in figure 4 (black lines and right axis) the 
Tln(t/ro) plot corresponding to the time magnetic relaxations of the ferritin sample. The 
scaling is reached using tq = 10 -8 s. One can see that the relaxation curves depart from 
the master curve at temperatures below 5 K, indicating the presence of quantum tunneling 
as the temperature arrives near the crossover temperature. Below ~2 K all the curves are 
parallel, showing the temperature independence of the quantum relaxation. 

We propose a new scaling law, equivalent to the Tln(t/ro) plot in time magnetic relax- 
ations, for the total change of magnetization of a volume distributed sample in a Landau- 
Zener process using a T\ia(t e ff /r ) plot, where t e ff = AH e ff/a. In figure 4 (open circles 
and left axis) is shown the Tha.(t e ff /r ) plot of -P^m = 1 — Pam, m order to compare with 
the magnetic relaxation master curve. The scaling is obtained with r = 10~ 8 s and AH e ff 
= 5 Oe. It is observed that the data collapse into a master curve for temperatures higher 
than ~5 K. The value of the effective resonance width, AH e ff = 5 Oe, is two orders of 
magnitude smaller than the width of the zero resonance observed in the magnetic hysteresis 
loops at the same temperatures, AH ~ 1000 Oe 0,|§-§[ and associated to thermally assisted 
resonant quantum tunneling P,B|,P| • The same phenomena was previously observed in molec- 



ular clusters |12|,|13| . In principle, the width of this resonance is associated to the quantum 
splitting of the blocking level, rag. This is the level through which the quantum tunneling 
occurs at a given temperature. In ferritin the width of the resonance is associated to the 
distribution of quantum splittings of the blocking levels due to the different volumes of the 
particles of the sample. This fact, together with the random orientation of the anisotropy 
axes respect to the applied magnetic field, makes the width of the zero field resonance to 
be several orders of magnitude higher than the width of the quantum splitting of one of 
the particles of the sample. However, the scaling law proposed here takes into account the 
effect of an average particle of the sample. Due to this, the physical meaning of AH e ff 
extracted from the master curve can be attributed to the width of the zero field resonance 
for an average particle of the sample. That is, we may associate AH e ff with the quantum 



splitting of the effective blocking level, A e //, of the distribution of particles in ferritin in 
the following manner: A e ff ~ gfigSAHeff. Using S ~ 50, we obtain A e jj ~ 700 MHz. 
Taking into account the uncertainties associated to the random orientation of the anisotropy 
axes of the particles it seems clear that the obtained value of the quantum splitting of the 
effective blocking level agrees with the ~1 MHz resonance found by Awschalom et al. |TJ 
and attributed to the quantum splitting of the ground state of ferritin particles. 

In conclusion, we have obtained temperature-independent magnetic relaxation from both 
magnetic viscosity measurements and Landau-Zener viscosity down to 100 mK. We have 
proposed a new scaling law for the probability to change the magnetization in a Landau- 
Zener process in a sample with barriers distribution. The excellent agreement between 
the two studied methods (time magnetic relaxation and Landau-Zener process) permits to 
establish that quantum tunneling is the process governing magnetic relaxation in ferritin at 
low temperatures. From the comparison of these two methods one can extract additional 
information about the magnitude of the effective quantum splitting playing the main role 
in the low-temperature magnetic relaxation of a sample distributed in volume. 

The authors thank the financial support from IST-1999-29110 and MEC Grant number 
PB96-0169. 
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FIGURE CAPTIONS 



Figure 1: Magnetic viscosity as a function of temperature of ferritin sample extracted 
from time magnetic relaxations at zero field. The inset shows the distribution of volumes of 
the sample (extracted from 

Figure 2: ZFC-FC magnetization versus temperature in ferritin recorded at H = 200 
Oe. The inset amplifies the low temperature zone. 1/T quantum superparamagnetism be- 
haviour is observed below 0.5 K. 

Figure 3: Probability to change the magnetization as the zero field resonance is crossed 
at different sweeping rates of the applied magnetic field and at different temperatures. The 
dependence with 1/a is clearly logarithmic. 

Figure 4: Tln(t e ^j/r ) plot for both time magnetic relaxations, t e ff = t (black lines, 
right axis), and Landau-Zener relaxations, t e ff = AH e ff/a (open circles, left axis). The 
values of r and AH e ff used to obtain the scaling are 10 -8 s and 5 Oe, respectively. The 
inset shows the temperature dependence of the Landau-Zener viscosity. 
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Figure 2: Duran et al. 
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Figure 3: Duran et al. 
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